On the torsion effects of a relativistic spin fluid in early cosmology 

G. de Berredo-Peixoto ^ * and E.A. de Freitas ^'^ ^ 

1. Departamento de Fisica, ICE, Universidade Federal de Juiz de Fora 

Campus Universitario - Juiz de Fora, MG Brazil 36036-330 

2. Colegio Tecnico Universitario, Universidade Federal de Juiz de Fora 
Av Bernardo Mascarenhas, 1283 - Juiz de Fora, MG Brazil 36080-001 



Abstract. In this work we investigate the effects of torsion in the framework of 
Einstein-Cartan theory in early cosmology. We study solutions for a homogeneous and 
isotropic relativistic Weyssenhoff spin fluid with dynamical timelike axial current, also 
homogeneous and isotropic. The general solutions can mostly be described by means 
of three particular solutions. The properties of these solutions (such as singularity 
avoidance and primordial or late accelerated expansion) are analysed and depend on 
the relations between the source parameters. 

PACS: 04.20.-q 04.50.Kd 98.80.-k 

Keywords: Torsion, Weyssenhoff fluid. Cosmology, Inflation. 



* Email: gbpeixoto@hotmail.com 
^Email: emanuel@fisica.ufjf.br 



1 



1 Introduction 



The effects of non-Riemannian structures in the framework of Cosmology have been studied long 
ago. The most simple generalization of General Relativity (preserving metricity) is achieved by the 
introduction of an asymetric connection, with torsion as its antisymmetric part. The Riemannian 
case is obtained by imposing zero torsion. It is possible to consider cosmological models in the 
framework of more general non-Riemannian structures, such as non-metricity (see for example the 
review by Puetzfeld [1]), but these cases will not be treated here. There are a large number of 
papers dealing with torsion in several approaches, with wide applications. For the introduction of 
the foundations of the theory of gravity with torsion, see Ref. [2], and for a more recent review, 
including the quantum aspects of torsion, see Ref. [3] . 

In the seventies, it was discovered that the singularity avoidance and inflation can be induced 
by torsion in the Einstein-Cartan theory, by Kopczynski [4], Trautman [5] and Hehl et al [6]. In 
the Einstein-Cartan theory, torsion is not dynamical and is completely expressed in terms of the 
spin sources [2]. Thus, in order to study the effects of torsion in the Einstein-Cartan theory, one 
has to introduce matter with spin. This can be done in several ways. One of the most natural ways 
is to consider, besides the Einstein-Cartan action, the Lagrangian describing spin-1/2 Dirac fields, 
minimally interacting with torsion. It is possible to describe this theory as a modified General 
Relativity with a spin-spin contact interaction. According to Ref. [7], torsion does not prevent the 
initial singularity, but rather enhances it. In addition, torsion can provide accelerated expansion 
phase in a metric-torsion theory with matter described by Dirac and Rarita-Schwinger spinors 
[8]. In this scenario, the authors of Ref. [9] described the transition between early accelerated 
expansion and decelerated one in terms of massive Dirac fields (see also, for example, Ref. [10]). 
It is possible of course to consider another Lagrangian for gravity and torsion, in a more general 
context of Poincare Gauge Theory of Gravity, with quadratic terms in curvature strcnghs. Within 
this possibility, one can mention for example the works in Refs. [11, 12] (see also references therein), 
which take into account a spinless matter. 

Another way to introduce torsion is to consider a fluid with intrinsic spin density, which in 
principle does not admit a Lagrangian description based on spinor fields. One has to postulate a 
spin correction to the energy-momentum tensor. For exemple, Szydlowski and Krawiec [13] have 
studied the cosmological effects of an exotic perfect fluid^ known as the Weyssenhoff fluid [15], as 
well as the constraints from supernovae la type observations, concluding that the dust Weyssenhoff 
fluid provides accelerated expansion but it can not serve as an alternative to Dark Energy. Also, 
in Ref. [16], Puetzfeld and Chen derived some constraints from supernovae la data in a different 
scenario of non-Riemannian geometry. Gasperini [17] considered the Weyssenhoff fluid with its 
energy momentum tensor (derivable from a Lagrangian) previoulsy improved by Ray and Smalley 
[18], with spin as a thermodynamical variable. In Ref. [17], torsion provides singularity avoidance 
and early accelerated expansion, but the expansion factor of the cosmological scale, a{t), is too 
small, unless the state equation parameter w {p = wp) of the spin fluid is fine tunned in a very 

^This fluid presents an intrinsic spin density (See also Ref. [14] where a quantum treatment was applied for this 
fluid in early cosmology) . 
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special way. See also Ref. [19] where similar analysis was performed in a more general context 
without assuming a particular metric. Obukhov and Korotky [20] formulated a more general 
variational theory describing the Weyssenhoff fluid and also applied to cosmological models with 
rotation, shear and expansion. Recently, Bohmer and Burnett [21] introduced a special spinorial 
matter satisfying the Cosmological Principle [22] (see also [23]), which, together with interaction 
terms and without cosmological constant, mimics vacuum energy responsible for inflation. 

It is worth mentioning that many works without torsion but with different matter sources 
(besides the scalar field) inducing inflation can be found in literature. For example, Golounev, 
Mukhanov and Vanchurin [24] proposed a scenario with massive non-minimally coupled vector 
field, which induces inflation. In homogeneous and isotropic universe, these vector fields behave as 
a minimally coupled massive field, and they can be introduced as an orthogonal triplet or in a large 
number randomly oriented in order to provide isotropic inflation. There are no clear understanding 
about the physical motivations for these vector fields. The authors of Ref. [25] consider also 
time-like non-minimally coupled vector fields (see also [26] where cosmological perturbations were 
studied) . Among theories without torsion, one can mention also Ref. [27] , where the authors show 
that the dark energy can be described by means of the usual electrodynamics with a non-linear 
adding term. Also, one can consider a timelike vector field, responsible for violation of Lorentz 
symmetry [28] . In the present work, the timelike axial current not only violates Lorentz symmetry, 
but it is related to torsion in the context of the Einstein-Cartan theory, and it is originated by 
the Weyssenhoff fluid. On the other side, wc can cite Ref. [29], where the possibility of inflation 
induced by non-standard spinors was investigated, as well as their imprints on CMB anisotropics. 

We consider the Einstein-Cartan theory with both the spin-spin contact interaction and the 
Weyssenhoff fluid. In practice, one has an axial current and a spin density as additional sources 
besides the usual perfect fluid energy density and pressure. In early universe, any matter content has 
very high temperature, so we let w fixed corresponding to a radiation fluid {w = 1/3), which mimics 
the ultra-relativistic regime. If the model is supposed to describe early inflation, then the spin fluid 
is an exotic form of matter which can play the role of vacuum energy. The matter content of the 
model is given by this exotic fluid (with w = 1/3) plus a timclikc axial current, homogeneous and 
isotropic. The anisotropic case is very interesting (providing description of primordial anisotropy), 
although it is beyond the framework of the present work. 

This paper is organized as follows. In section 2, we shall present a brief introduction to Einstein- 
Cartan theory, including the dynamical equations of the model. In section 3, we investigate the 
solutions of the model for different cases. In section 4, we draw our conclusions and final remarks. 

2 Einstein-Cartan Theory and dynamical equations 

The action in the Einstein-Cartan framework is given by 




(1) 
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where metric has signature (H ), = IGvrG (we use units such that h = c = I), A is the 

cosmological constant and R is the Ricci scalar constructed with the asymmetric connection^, T'^a/j, 
which, by using the metricity condition iy adixu = 0) and the following definition of torsion, 

can be expressed as 

f'^a/3 = r^a/3+i^^/3, (2) 

where T^^ctiS is the Riemannian connection (Levi-Civita connection) and the quantity K^^^p is the 
contortion tensor, given by 

The term Cm is the Lagrangian describing matter distribution. We consider here the following 
matter Lagrangian: 

= ^AC + ^SF , (3) 

where Csf is the Lagrangian of the spin fluid [18] and Cac is the external source, present in the 
minimally coupling Dirac sector (see, e.g., [3]): 

Cac = J''S^. (4) 

Here, 5^ is the axial part of torsion, defined by = exp^^T^P" {^Xpcru is the Levi-Civita tensor, 
with £oi23 = sZ—g), and is the external axial current^, J'^ = < 'ipj^j'^ij; >, where this average 
is due to quantum effects (see [3]), such that is a vacuum property, responsible for Lorentz 
violation (see also [30]). 

In order to vary the action and get the dynamical equations, one has to define what are the 
independent variables. We choose g^'^ and T"^^ as independent dynamical variables, and as 
a quantity defined by the symmetry violation of the vacuum, which is, as wc shall see, dependent 
from dynamical variable g'^'^. The spacetime metric is the spatially flat homogeneous and isotropic 
metric such that 

ds'^ = dt'^ -a{t)'^{dx'^ + dy'^ + dz'^). (5) 

It is natural to assume that, in comoving frame, is a homogeneous and isotropic vector, otherwise 
it would break the isotropy of the universe. Thus, is a timelike vector so that we let J'^J^ = J^{t). 

^AU quantities with an upper tilde are constructed with the asymmetric connection, and the corresponding quan- 
tities without tilde are constructed with the Riemannian (symmetric) conection. 
^The matrix 7^ is the chiral Dirac matrix 7® = (i/4!)e"'^'"'7a7/37^7,, = 170717273- 
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2.1 Variational principle and dynamical equations 

The dynamical equations for metric fields and torsion, in terms of the sources, can be obtained 
respectively by the usual procedure'^ 

- ^ ^^<^ . (a) 



vo. 



where 



and 



^ ._ 2 (5(v^(£aC + ^5f)) ._rrAC , rrSF 
Q va 

1 5{^{Cac + Csf)) ,^ 

\ va , ( Q \ va 

C'm — -7= •= {f^ACju. + [(^SF)ix ■ 

V g va 

It should be noticed that instead of procedure (7), one can verify that 

^ = gives 1 = (8) 



where 



1 5{-^/^Cm) 

g "^ixva 



fj-ua . fiva 
T An ~r T, 



Equation (8) is totally equivalent to equation (7) , thus we shall use it for the convenient correspon- 
dence with notations in literature. 

By the convention for the curvature tensor in the form'' R^xajS = 2r'^x[p,a] + ^''xpi^'^pa — 
Vf'xaT^pjB, one can achieve, disregarding total derivatives, the relation 

j d^x^R = J d^x^ (r + K'%^K''\ - K'^xK''^) 
Variation with respect to T^^a/S gives 

rp^a^ + 25^["r^] = kV^"'' , (9) 

where = T^^p. In order to express r^"'^, let us mention that the contribution from spin fluid is 
given by Tgp^ = —^S^'^u^ [18], where S^"' is the spin tensor (antissymmetric) , and is the fluid 
four- velocity^. By straightforward algebra, we can obtain the expression for T^"^^ = + T^p^ 
and consequently. 



(10) 



^We use similar notations from literature, e.g., Rcf. [17]. 

''We use the symbols [] and () to denote antisyinmetrization and syrnrnetrization, according to 

-^li^i^] = — A^ij) and -^{^iv) = gC^M"' + ^■'m) ■ 

^In previous papers, like, e.g., Ref. [17], the expression Tgp^ has the opposite sign. The reason is that we adopt 
LsF with different sign in order to reproduce the same fluid dynamical equations usual in literature. 
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Using the Weyssenhoff condition^, S^'^Ua = 0, from the above equation we can derive 

Tf^f^ = -K^ !^2e'"^'^PJp - ^S'^'^u'^^ (11) 

and 

Sa = Wa(2if'^"^) = 12acV, + ^K%apaS'"'u'^ ■ (12) 

Variation with respect to g^^'^ can be done in a straightforward way, giving 

ixv yfxu^^ ^ fj,a^ vX 2 Pl^ ^'^ 2 1/ ixpX-'-v 

+ ^g^, (AK",aK^''x + 2T"''%xa - T^^^T^^a) = yT^. • (13) 

To write the dynamical equations, one can substitute (11) and (12) into the appropriate quantities 
in (13), including T^i,. In this way we can rewrite equation (13) in the form 



4(2 1 A 1 A 



SF\ 
liv I 1 



(14) 



where = J^J'^- Notice that from (11) one has T"^^ = 0, such that several terms in (13) vanish. 
Thus, with axial current and the spin fluid satisfying the Frenkel condition, S^'^Ua = 0, there are 
only tracclcss degrees of freedom of torsion. 

The next step is averaging the above equation. For this purpose, a natural and simple assump- 
tion is < Saf3 >= 0, which means that although the spin tensor might have a particular direction 
in the microscopic scale, its mean value vanishes at macroscopic domain (i.e., the particles have a 
randomic spin distribution). Let us define (see Ref. [17]) 

< Saf^S'^f^ >= 2a^ , (15) 

such that 

2 

< Sf^^SuX >= i^igfii^ - u^u^)a^ . 

Now we have to express T^^ and T^^ in terms of the sources. For T^jf , one can obtain the 
formula (see, e.g., Ref. [17]): 



1. 



+ -T^^^^Sf" p^u^ + {(p + p)uaup - pg^p} . (16) 
Substituting (11) into (16), one obtains 



2 2 

+ 2k^ {{p + p)Uo,Uj3-pgal3} ■ (17) 



^Also known as the Frenkel condition. It is included by hand, but emerges automatically in the formalism proposed 
by Obukhov and Korotky [20] . 
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We achieve, by averaging, 

- -^UaUpa^ + —gapCf'^ + {{p + p)UaUfi - pQap] ■ (18) 

It should be mentioned that, although < S^y > vanishes, it is possible to take < Sfj,uJ°' >^ 0. 
Nevertheless, we argue that the assumption < Sf^i.J" >= is correct because S^i^ is randomic only 
in the 3-space, and has not any spatial component. 

Let us finally express T^'^ in terms of the sources. The variation of ^/—gCAC with respect to 
g^^'^ should be done with special care, since (and also J^) depends on g^'^. This computation is 
done in the Appendix. The result can be written as 

T^i^ = -2£(^°'^T^)c«/3'^p - S(nJ„) . (19) 
Substituting the sources (11) and (12), we have 

= -Sac'^m- J' - 4acV^ J, - Kh^^^P\,)SfSpJx + y e(/''V,)5pAn;3 . (20) 
Now, one can obtain the following dynamical equation after averaging (using < J" >= 0): 

f 1 1 1 

Gij,u = ti^ { -'ig^.vJ'^ + -^9tivcr^ - -^u^Uya"^ ^ + y {(p + p)y'tiUy - pg^,v} + Ay^v , (21) 

Dynamical equations 

Let us consider the fluid in a relativistic regime, such that p = p/3. For the metric (5), the relevant 
components of equation (21) can be written as 

where the dot above variables means time derivative. From these equations, one can get 



It is remarkable that the axial current enters in the above equation with different sign from the spin 
contribution. In this equation, the axial current acts as a kind of density energy of some ordinary 
matter. We see that only the cosmological constant and the spin contribute to the accelerated 
expansion. 

The energy conservation law can be obtained by comparing equation (23) with the time deriva- 
tive of equation (22): 
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Let us remark that the above equation describes energy conservation in the general case when all 
fields interact with each other. Of course, the spin contribution (a^) always interacts with the fluid 
itself, because spin is a fluid property. On the other hand, we have freedom to choose an interacting 
axial current and a non-interacting one. These possibilities define two classes of solutions. The 
non- interacting axial current case was studied previoulsy in Ref. [31]. The case without spinning 
fluid can be found in Ref. [3], where a different variational procedure was realized. 

Equations (22) and (25), along with the initial conditions p{to) = po and a{to) = oq = 1 {to is 
the present time), determine the dynamical solution of the model. For the CMB, we know that 

po ~ Pc^rad- As 

o rr2 • 

f^rad ~ 10^^ and pc = - 4 X W-^\GeV)'^ (with H = -) , 

one can estimate po ~ 10~^^(GeF)^ for pure radiation. Of course, this value does not necessarily 
correspond to the exotic spin fluid, but can serve as a reference (or upper bound). In what follows, 
we shall consider in most cases a much lower density for the spinning fluid, say, po ~ 10~^^(GeF)'^. 



3 Solutions with interacting axial current 

In searching for solutions, we should specify the dependence of and on the density, p. One 
can adopt = 7/9^^^ (7 = positive constant) [32]. We identify the axial current as coming from 
the spinning fluid. In doing so, the natural choice for is = /3p^/^ (/? = positive constant). 
With these assumptions, the dynamical equations for the model are 

12K2d7p3/2 - Uap ^26) 



16a - 3K2a7/9V2 _ lUH?aPp^/'^ 
and 



P+|-kV3/2(^/3 + ^7). (27) 



'k2 a 



The above system is very complicated and can not be solved analytically. However, we can 
extract relevant information just from (26). Let us rewrite it: 

dp ^ 12/^2^^3/2 _ 

da 16a - 3k2o7pV2 _ UAH^af^p^l"^ ' ^ ' 

For a particular p = pf, there is a fixed point, i.e., dp/da = 0. When p(a) reaches pf, it ceases to 
vary. One can express pf as 

_ 256 

Also, there is an apparent singularity when p = pc such that dp/ da — ^ 00. The value for pc is 
obtained by straighforward algebra: 

_ 256 

~ (3^27 + 144«;2/3)2 • 
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Actually, the system is not only apparently singular at pc, but it is really singular, because a (from 
(27)) can not vanish for pc. Notice that always pc < Pf- Let us investigate the three unique 
possibilities: (i) po < pc < pf, (ii) Pc < Po ^ Pf and (iii) Pc < Pf ^ Po- For the first possibility, 
we have p' < in the whole interval p < pc- Thus, as time goes backward, density is increased 
until it reaches pc- It means that in some finite a > in the past, p' is infinite. As it has not clear 
meaning, we must reject the possibility (i). 

For the second possibility, we have p' > in the whole interval pc < p < Pf- Using similar 
reasoning as before, p' becomes again infinite for some finite a in the past. In the last possibility, 
p' < in the interval Pf < p, such that p never reaches pc- Thus, the range given under (iii) 
appears to be the only viable physical choice. 

By the above considerations, the following inequality must be satisfied: 

256 

In numbers, it means that the dimensionless parameter 7 has a lower bound determined by pQ. 
For Po = 10-5^ GeV^, we achieve a very large lower bound^: 7 > 1.58 x 10^^! This lower bound 
imposes the lower bound cr^(io) = 7Po^^ > 1.6 x 10"-'^''' GeV^. 

However, a carefuU analysis shows that the assumption cr^ , oc p^/^ is not rigorous. Indeed, 
one should start from a more simple and fundamental assumption, by considering the dependence 
of cr^ and (or {i^ipY) on the scale factor, as 

(7^ oc a-^ . (29) 

Notice that in the previous paper [17] the fluid density satisfies p oc a~^, such that cj^ oc p^/^. 
Nevertheless, p oc a^*^ is clearly not an exact solution of the equation (25). Thus, the ansatz (29) 
will be adopted from now on. Let then 

J2 = 4 and (t2 = 4- (30) 
The equations of motion (22) and (25) can be written in the form 

a^- " 48a6;+ 6^+3' ^^^^ 

dp ^ _4p _ 36^ ^^2) 

da a a7 

It is quite difficult to solve analytically the above system, so it is convenient to treat it numerically. 
The first observation is that the quantity dp /da in equation (32) is negative and lower than —Ap/a, 
thus the effect of the axial torsion is that, as far as po (= p{to)) is taken to be the same quantity 
for both cases (with and without axial current), the values for p are higher for t < to- 

Now, let us consider equation (31). Its right hand side must be positive for all values of a{t). 
Choosing the special case a = 1, we arrive at the following condition (for po = 10~^^ GeV^): 



^Notice that A = 5 x 10"^* GeV^ and = 3.38 x lO'^^ GeV^ 
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which defines an upper bound constraint for the source parameters Jq and (Tq. In equation (31), 
the term proportional to has the same dependence on the scale factor than the non-Riemannian 
parameters in Ref. [16]. Although a different model is studied in this work, this coincidence opens 
the possibility for estimating some constraints from supernovae la data applied to the present 

2 

model^. As a result, we get Jq + Is — ^'^^ ^ 10~^^GeV^, which is remarkably similar with (33). 

Before investigating general solutions for (31) and (32), let us consider the obvious particular 
solution of (32): 

p{a) = g . (34) 

Substituting this solution into (32), one achieves the constraint 

po = 18k . (35) 

Notice that in this case the equation (33) establishes effectively an upper bound for ctq , since Jq is 
attached to po- The equation (31) reads 

= ^ + A „,e.e .:=„'(24-|). (36) 

There are three solutions depending on the sign of 6: 

( [39 V^' 
a{t) = S Y ^ sinh {VSAt) \ for > , (37) 

where the integration constant was already fixed by a(0) = 0, and the second solution is given in 
the implicit form: 



a(i)3 + y^ + a(i)6 = y-^exp(V3Ai) for e<0. (38) 

In the above solution, the integration constant was choosen such that a(0) = amin = (— 3^/A)-^/^, 
with amin being the minimum value of a which can be found by the condition 6 + Aa^ > 0. Thus, 
for ^ < 0, there is a singularity avoidance, and the universe undergoes an accelerated expansion all 
the time. In contrast, solution (37) does not prevent the model from having initial singularity. 

The third solution comes from the case 6 = 0. It is given by a{t) oc eV^^ (De Sitter). It 
is remarkable that in this case the Universe expands as it was empty with only the cosmological 
constant, but there is a fiuid with p oc a~^. 

It is interesting to extract information about the source parameters, Jq and ag, just from 
the experimental constraints such as the known age of the universe^^, to = 13.7 billion years 
= 6.56 X 10^^ GeV~^. Let us consider, for example, the solution (37). By using the expression for 
a{t) in (37) at the equality a{to) = 1, one can get 6 = 4.18 x 10~^^ GeV^ (similar quantity can be 
found for \9\ in the case (38)). Now, using equation (24) in the form 

d__29 A 
a 3 ' 



^We thank the anonymous referee for this indication 
"'Notice that this proceedure is not rigorous if we rem 
by the conventional matter content, such as perfect fluid in dust form 



^"Notice that this proceedure is not rigorous if we remember that in late times, the role of torsion will be supressed 
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we achieve, by direct substitution, 

o 3 I J 

It is possible to show that this quantity is positive if t > 2.96 x 10^^ GeV~^, what means that, 
taking (9 = 4.18 x 10"^^ GeV^, the expansion of the universe is accelerated for a > 0.61. 

Notice that in the absence of the cosmological constant, the quantity 9 must be positive, admit- 
ting then the solution (37). Only a non-null cosmological constant can provide singularity avoidance 
(for Jg 7^ 0). It is remarkable that the theory with axial current and spin fluid, satisfying (35), 
requires the introduction of A, otherwise the model would have no solutions for all possible values 
of CTg. For Jg = 0, the model reduces to the one studied by Gasperini [17], with early and late 
accelerated expansion, and also singularity avoidance. In this case, the early accelerated expansion 
takes place in a very short period of time, and it should be mentioned that Jq 7^ can be taken 
as an important generalization, in a Lorentz violating theories or as a vacuum quantum effect, for 
example. 

In fact, condition (35) seems to be quite particular. However, as we shall see below, the general 
solutions for /3(io) / 18k^ Jg can be described by the particular solutions dictated hy p oc a~^. 

General Solutions 

According to previous considerations, we know that p oc is a particular solution, which demands 
(35). Indeed, one can perform numerical integration of equation (32), using the Mathematica 
software, starting from the point p{aQ = 1) = po satisfying (35). The integrated curve for p{a) will 
be exactly the curve for p oc a~^, which can be drawn as a straigh line with negative slope (—6) in 
the logarithm scaling. 




a 



Figure 1: The integral curves, in logarithm scale, for p(a), from initial conditions p{l) = 10 
GeV^ and p{l) = 18k = 6.08 x lO"^'^ GeV coincide, except in the smah region shown in the 
right plot. 

In order to numerically integrate, we choose some Jg compatible with condition (33), say, 
Jg = 10~^^ GeV^, and po to be different from (35). Let us choose a lower value: po = 10~^'^ GeV^. 
The integration is shown in the left side of Figure 1, together with the integration obeying (35), 
corresponding to /) oc a~^. It is remarkable that both curves coincide, but it is essential to stress 
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that they coincide in a certain (and wide) interval. The same plotting in the right side of Figure 1 
is drawn in an interval much more closer to a = 1, where the difference between the two solutions 
becomes clear. 

It seems that this behavior is universal: the curve p{a), for any po, differs from the solution 
p oc in some region close to a = 1, but in the remaining region, say, < a < 0.1, the solution 
is very close to p(a) = ISk^Jq/q^. Now, remind that there is no much room for big values of po, 
and we are mostly considering early cosmology, thus the particular solutions discussed previously 

are quite general and instructive. 

We should mention that the case Jq = [17] is completely different, because p(a) will be 
substantially affected, such that p oc a~^. 

4 Conclusions and Discussions 

We have investigated the cosmological effects of the axial current together with the relativistic 
spin fluid (p = p/3) in Einstein-Cartan theory. As already known from literature [4, 5, 6], torsion 
provides singularity avoidance and accelerated expansion. The contribution from axial current, 
however, favours a decelerated expansion, in contrast to the spin fluid. 

There are two classes of solutions: one with an external non-interacting axial current, and other 
with interacting axial current (i.e., time-dependent axial current). The first one was considered in 
previous works (see [3] for the case without spin fluid and with a conformal global axial vector, and 
[31] with a global constant vacuum axial vector and a spin fluid). 

In the present work, the axial current is assumed as a composite field, J'^ =< 'ip'y^^'^'ip >, where 
the Dirac fields (presumably) describe the fluid itself. The axial current interacts with the spinning 
fluid, as realized by the energy-momentum conservation. It is essential that this feature determines 
the dependence of energy density on the scale factor, a, which is substantially different from the 
case with non-interacting axial current. It is natural to assume that both spin and axial current 
decay in similar way as the universe expands. 

We conclude that, based on the ansatz oc and oc a~^, the general solutions behave 
(in a relevant domain) as the particular solutions coming from p oc a^^, such that their properties 
are the same. These properties were analiscd in details for the particular solutions: If the source 
parameter 6 is positive, the universe has an initial singularity and its expansion is decelerated 
(a < 0) until some epoch, when late accelerated expansion begins. This epoch depends basically 
on the parameter 6, which can be determined from the known age of universe. By doing this, 
late accelerated expansion starts at a = 0.61. In both particular solutions, Jg can be fixed by the 
present value of density, po. Thus, ctq determines the sign of 9. 

For the case 9 < 0, universe has an accelerated expansion all the time, and the solution features 
singularity avoidance. In this case, equation (33) represents an upper bound for the quantities Jg 
and (Tq. 

As discussed above, the general solutions have a remarkable behavior (as the shift shown in the 
right side of Figure 1) in the region close to the present day. This is very strange in the physical 
point of view, because it would be an enormous coincidence if the shift of the actual curve occurs 
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right on the present time, to. Thus, it seems that the physically reasonable solution must be the 
particular one, p oc a~^. As a consequence, the parameter Jq should be related to the present 
density by equation (35). 
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Appendix: veiriational proceedure 

Here we shall calculate the functional derivative of J'*. To do so, one must take into account that 
= -^-yS^M^^ with 7^ = (^/4!)£°^'*''7Q,7^7^7^^. As 7^7^* = 4, we can write 5{'^al°') / ^9^^ = 0. With 
this equation, one can express Sja/Sg^'^ in terms of Sj/^/Sg^^ and vice- versa. Using = g''^^x, 
we achieve 

SjP 1 (57 1 

S^ = 2^i,M ^ = --5^(^7,). (39) 

Now, we know that £°'I^p^ = E'^/^P^ /^/^, where E°'l^''^ is the Levi-Civita symbol, which is 
clearly independent on the metric. Thus, we obtain 

-3^ = 2^^'^^ -^ = -^g,.s^ppx. (40) 



Using (39) and (40), one can get 



= 0. 



SgP'" 

With all these results, one obtains by straighforward algebra the variation of J" and Jq.: 

Similar computations can be performed for = e'^^P^T^pp = e'^^P^gaaT" pp. After all, the 
variation of \/^J^Sfj^ can be expressed by means of result (19). 
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